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Abstract 

We describe a family c f>\ of dynamical systems on the unit interval 
which preserve Bernoulli convolutions. We show that if there are pa¬ 
rameter ranges for which these systems are piecewise convex, then 
the corresponding Bernoulli convolution will be absolutely continuous 
with bounded density. We study the systems (j)\ and give some nu¬ 
merical evidence to suggest values of A for which cf>\ may be piecewise 
convex. 


1 Introduction 

In the study of self similar measures corresponding to non-overlapping iter¬ 
ated function systems, there is a natural way of defining an expanding dy¬ 
namical system which preserves the measure and which allows one to study 
various properties such as dimension. The case of self similar measures with 
overlaps is much more involved, and it is not clear how best to study them 
using dynamical systems. 

Bernoulli convolutions are a particularly well studied family of self-similar 
measures. For each A G (0,1) we define the corresponding Bernoulli convo¬ 
lution v\ to be the distribution of the series 

OO 

(A -1 — 1) ^ a*A* 

i =1 


l 



where the digits a* are picked independently from digit set {0,1} with prob¬ 
ability Equivalently, Bernoulli convolutions are the unique probability 
measures satisfying the self similarity relation 


= ^(>a ° T 0 + u x o Ti), 

where the maps Tj : M. —>• M. are dehned by T,(x) = | — (A -1 — l)i. For 
A G (0, |), the self similar measures are generated by a non-overlapping 
iterated function system, and are invariant under the interval maps (j)\ given 
by 

( A _1 x x G [0, A] 

4>\(x) = < 1 X G (A, 1 - A) . 

[ 1 — A -1 (x — (1 — A)) x G (1 — A, 1) 



Figure 1: The maps 4>x f° r A equal to 0.2, 0.4 and 0.5 respectively 

The main aim of this article is to extend the definition of <p\ to the overlapping 
case, when A G (|, 1), and to study v\ using these interval maps. 

There are a number of long standing open questions about Bernoulli con¬ 
volutions, chief among which is the question of for which parameters A the 
corresponding measure v\ is absolutely continuous. It is known that each 
Bernoulli convolution is either purely singular or absolutely continuous, see 
[3] . If A is the inverse of a Pisot number then v\ is singular, see [5] , and in fact 
has Hausdorff dimension less than one, |9j. In [HJ Garsia gave a small, explic¬ 
itly dehned class of algebraic integers for which v\ is known to be absolutely 
continuous. Solomyak proved in |13j that v\ is absolutely continuous for al¬ 
most every A G (|, 1), and in [12] Shmerkin proved that the set of A G (|, 1) 
admitting singular Bernoulli convolutions has Hausdorff dimension 0, but the 
question of determining the parameters A which admit absolutely continuous 
Bernoulli convolutions remains open. 
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There are other interesting open questions regarding Bernoulli convolutions. 
For example, is it the case that any singular Bernoulli convolution must have 
Hausdorff dimension less than one? Do there exist intervals in the parameter 
space for which every Bernoulli convolution is absolutely continuous (and 
even has continuous density)? Does the density evolve continuously with A? 

Similar questions exist in the study of invariant measures associated to var¬ 
ious one parameter families of interval maps, and in this area a good deal 
of progress has been made mu ITU] , With this in mind, we extend the 
definition of the generalised tent maps (p\ to the overlapping case. These 
tent maps preserve the corresponding Bernoulli convolutions v\. They are 
described implicitly in terms of the distribution F\ of v\, and while we are 
able to write down explicit formulae for the (p\ only in some special cases, 
we are able to prove some general properties. 

In particular, we prove that if (j)\ is piecewise convex for all A in some interval 
(a, b ) then the corresponding Bernoulli convolution is absolutely continuous 
with bounded density. For each x G [0,1] the map x t->- <j>\{x) is continuous 
in A, and convexity is preserved by passing to limits in a continuous family 
of functions. Thus, piecewise convexity of the functions (J)\ seems like an 
appropriate vehicle for passing from almost everywhere absolute continuity 
to everywhere absolute continuity for parameters in certain ranges. We can 
show that 0 a is piecewise convex for certain special cases, and remain opti¬ 
mistic that one may be able to prove analytically that the map <fr\ is piecewise 
convex in certain parameter ranges. For the moment however, our results 
on piecewise convexity are restricted to some special values of A, although 
we are able to run numerical approximations for any A. There have been 
previous numerical investigations into Bernoulli convolutions, we mention in 
particular the work of Benjamini and Solomyak [1] and of Calkin et al Eli- 

in the next section we define the maps (j)\ in which we are interested and 
prove that they preserve Bernoulli convolutions. We prove some elementary 
properties of the maps (f>\ and give the maps explicitly in some special cases. 
In section 3 we prove various properties of v\ that would follow from (j)\ 
being piecewise convex, and in section 4 we give some numerical evidence 
on the piecewise convexity of (p\. Finally in section 5 we state some further 
questions and conjectures. 
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2 Generalised Tent Maps 


Let F\ : [0,1] —> [0,1] be the distribution of u Xl i.e. F x (x) : = i/ A [0,x]. F\ 
is strictly increasing because u x is fully supported. We define a map <f> x : 
[0,1] ->[0,1] by 

( )= ( F~\2F x (x)) *6 [0,1] 

M ’ l F-\2F x (1-x)) X6[l,l] • 

Since F x is strictly increasing on [0,1], the map <f> x is well defined. We will 
see later that 0 A preserves u x . 





Figure 2: Graphs of 0 A for A = 0.6,0.7 and 0.8. 

The map (j) X will be the chief object of study for this article. Since F x can 
be well approximated numerically, to known levels of accuracy, one can gain 
good numerical approximations to the maps (j) X . Three such approximations 
are displayed for different values of A in Figure [2j 

We begin by observing some simple properties of (j) X . 

Lemma 2.1. The map (ft x has the following properties for A G (|, 1). 

1 . 0 A (O) = 0 A (1) = 0. 

2 . *,©=!■ 

3. 4> x is strictly increasing on [0, |] and strictly decreasing on [1, l]. 

4■ <j>x(x) = (f x (1 - x). 

5. (p x is continuous. 
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Proof. We have that F x ( 0) = 0, Fa(|) = | and F x ( 1) = 1 because v\ is 
supported on [0,1] and symmetric about the point 1. Then points 1 and 2 
follow immediately. 

Part 4 can be seen to be true by looking at the piecewise definition of <f>\. 
Because v x [a, b] > 0 for each 0<a<6<lwe have that F\ is strictly increas¬ 
ing. Consequently (f> x is strictly increasing on [0, |] (and strictly decreasing 
on [|,1]). 

Finally, we observe that continuity of <f> x follows from the fact that v x is non- 
atomic and that v x [a, b] > 0 for each 0 < a < b < 1. Then both F\ and F^ 1 
are uniformly continuous, and so <f> x is continuous in x. □ 

We call maps satisfying the above properties generalised tent maps. Our first 
theorem is the following 

Theorem 2.1. Let A e (|, 1). Then v\ is invariant under (f>\. 

Proof. It is enough to show that for each a G [0,1] we have that 

z/ A [a, 1] = z'aO/’a 1 ^,!])- 
To prove this, we note that 

‘/’aMM] = [b, 1-6] 

where b G [0, |] satisfies 4>\{b) = a. But then 
u x [b,l-b\ = 2 v x [b,^\ 

= 2 (F x (±)-F x (b)) 

= 1-2 F x (b) 

= F x (l) - F x (Mb)) 

= F x (l)-F x (a) 

= v\[a,l], 

as required. □ 

Thus, if v x is absolutely continuous we see that <f x has an acip. We have not 
been able to prove the converse statement, that <f> x does not have an acip in 
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the case that v x is singular, this would be a useful statement which would 
make the relationship between the study of <f> x and the measures v x a little 
more straightforward. 

The following theorem shows that the maps <f x evolve continuously in A. 

Theorem 2.2. For each x G [0,1], A 0 G (|, 1) we have that <f x {x) —> 4> Xo (x) 
as A — y A 0 . 

Proof. Fix Ao G (|, 1). We rely on three facts for this proof. 

Firstly we use that the function F^ ] is continuous in x: for all 62 > 0 there 
exists 61 > 0 such that 

\x — y\ < 2ei =► \F~\x) - F i\y)\ < 62 . (1) 

Secondly we use that for each x G [0,1] the function F x (x) is continuous in 
A: for all e 3 > 0 there exists <5i > 0 such that 

|A - Ao| < S 1 => | F x (x) - F Xo (x) | < Cl . (2) 

Finally we use that for each x G [0,1] the function F A ” 1 (x) is continuous in 
A. For all e 3 > 0 there exists a d 2 > 0 such that 

|A - A„| < S 2 => | F?(x) - F£(x) I < £ 3 . (3) 

We £x x and let 5 = min{Ai, <5 2 } and |A — Aq| < 5. Then 


\Mx)-K(x)\ = K'(2F x (x)) - F^(2F^(x))\ 

< sup |F a '‘(2F a (i)) - F a _ 1 (j/) 

2F x (x)-2e 1 <y<2F x (x)+2e 1 

< \F-\2F x (x))-F- 1 (2F x (x))\+e 2 

< £3 + e 2 , 

Here the second line holds since 2F Xo (x) G (2 F x (x) — 2e,2 F x (x) + 2e) by 
equation [2] Then the third and fourth line follows from equations [l] and [3] 
respectively. Since e 2 ,63 were arbitrary, we are done. □ 

In the case that one knows the distribution F x , one can write down the map 
( t> x explicitly. In particular, for A = 2“^, it is not difficult to write down <f> x . 
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Example 2.1. In the case \ = ^, F\ is given by 




WH a + W-W 

l-(l + |V2)(l-i) ! 


Consequently <f>\ is given by 


x G 
x G 
x G 


0, 


1+72 


_%/2_l 

1+72’ 1+72 J 
I! 

1 + 72 ’ 4 


\[2x 


x G 

(1 + V2)X + 2_|_2v / 2 

x G 

1 

to 

^i/2-,y/ 2 

v 1+72 ) 

x G 


0, 


2 + 72 J 
1 1 


2+72’ 1+72 


1+72’ 2 


which is extended to the whole interval I\ using the symmetry around |. We 
have drawn the graphs of F\ and <f>\ in Figure [3| 


2.1 Further properties of cf>\ 

While we cannot write down explicitly, we can describe the behaviour near 
x = 0 and the rate of the blowup at x = The following lemma describes 
<f>\ near 0, and hence also the behaviour near 1. 

Lemma 2.2. We have that 


0 A (x) = A 


for x G [0,1 — A]. 
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Proof. Self similarity of the measures v\ give that 


F x (x) = ^ (F x (X-'x) + F x (X-'x - (A -1 - 1))) 

Then 

F x (fx(x)) = 2F x(*) = FxiX-'x) + FxiX-'x - (A" 1 - 1)). 
But because F x (x) = 0 for x < 0, we have 

F\(X~ l x — (A -1 — 1)) = 0 

for x < 1 — A. Then 

Fx(M*)) = 2F ^) = F x( ^x), 

for x G [0,1 — A], which completes the proof. 


( 4 ) 


It remains to find fx(x) for x G [l — A, |], and then by symmetry to define 
fx on [|, 1], We can also describe the nature of fx around x = \ for typical 
A. 

Lemma 2.3. We have that 

f 1 A _ log A 

f I- x ~ 1 — CX lo E 2 

for small x, where c is a constant that depend continuously on A. 




Conversely, equation [4] gives that for small 5 


W-S) = ^(1 + Fx(A _1 (l — 6) — (A -1 — 1))) 


giving 


1 - F a (1 - S) = 4(1 - F X ( 1 - A' 1 *)). 


( 6 ) 


Now suppose that <f>\ (| — e) = 1 — 5 for some fixed e and 5. Then by 
equations [5] and [6] we have that 



Iterating, we have 



and we see that we have a blow up of the form 



with c = fc lo § A / lo g 2 depending continuously on A. 


□ 


3 Piecewise Convexity 

Numerical approximations of the maps <f>\ suggest that there are ranges of A 
close to 1 in which the maps <f>\ are piecewise convex. For the value A = ^, 
one can see directly from the calculation of the previous section that <f>\ is 
piecewise convex. 

In this section we prove various properties of v\ that would follow from 
(f)\ | [q ; i] being convex. We use the term ‘piecewise convex’ as shorthand for 
the statement that (p\ is convex on each of the two intervals [0, |] and [|, 1]. 
The following theorem shows the relevance of the piecewise convexity of (f)\ 
to the study of Bernoulli convolutions. 

Theorem 3.1. Suppose that there exists an interval ( a,b ) C (|, 1) such that 
(p\ is piecewise convex for each A in ( a,b ). Then for each A G (a, b) the 
Bernoulli convolution v\ is absolutely continuous with bounded density. 
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We stress that if 0 a is piecewise convex for almost every A in ( a,b ), then it 
is piecewise convex for all A in (a, b), since the maps 0 a are continuous in A 
and convexity is preserved by passing to continuous limits. 

Proof. This theorem relies on results of Rychlik IIIl- 

Given a function g : [0,1] —>■ M, we define the total variation of g by 

n 

var g := sup Y 1 9( x i) - g( x i-i)\- 

0=XO<Xl<---<X n = l 

The function g is said to have bounded variation if varg < oo. Suppose that 
T : [0,1] —» [0,1] is a piecewise continuous map, such that there exists a 
function g of bounded variation satisfying g = 1 /1W' | almost everywhere. We 
consider the transfer operator L defined on functions of bounded variation 
by 

L f(x) = Y 9{y)f{y)- 

T(y)=x 

We put g n — g ■ (g ° T) ■ ■ ■ (g o T n_1 ). Then 

L n f{x) = Y 9n(y)f(y)- 

T n (y)=x 


Let C n denote the (n — l)th refinement of partition {[0, ^], (|, 1]} by T. 

In HU, Rychlik proved that 

var L n f < k var / + D\\f\\i, (7) 

where k = sup g n + max 6; var c „ g n and D = max Cn var Cn g n /\C n \. 

We can apply this to our tent maps, replacing T with <p\. Suppose that 
the tent map is convex on each of the intervals [0, |] and [|, 1]. Then 0 a is 
differentiable everywhere except for at most countably many points, and this 
derivative is increasing on [0, |) and on (^,1]. So there exists a function g 
which is of bounded variation, which satisfies the assumptions of HU, and 
which satisfies g = almost everywhere. We have 

sup g = g( 0) = g(l) = A and var^ij g = var ( i^ g < A, 
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with equality if and only if \(p' x (x)\ —> oo when x —> \ (which is the case by 


Lemma 2.3). From this we get that 


supg n = A n and var Cn g n < 2 n 1 X n , 


Combining this with ([7]) we get that 


var 


on-1 \n 

L n f < 2A n var / H- : — 

nun \ C n \ 


( 8 ) 


In the setting of our tent map, C n corresponds to the cylinders of generation 
n, and so C n depends continuously on n. In particular, for each n G N the 
value of min \ C n \ corresponding to (ft\ is continuous in A. 

By Rychlik [TT] , there is a unique non-negative function h\ of bounded vari¬ 
ation, such that \\hx\h = 1 and Lh\ = h\. The function h\ is the density of 
the unique absolutely continuous invariant measure of (ft\. If we pick n such 
that 2A n < 1, then (J8]) implies that 

2 n_1 A n 

var h\ < 2A n var h\ H- : ———, 

mm |C n | 


giving 


Hence we have that 


var h\ < 


2 n -'A n 1 
min C n 1 — 2A n 


sup h\< 1 + 


2 n_1 A n 1 
min \C n \ 1 — 2A n ’ 


and so h\ is bounded. Furthermore, since all of the quantities involved are 
continuous in A, there is a uniform bound on sup h\ across all of (a, b). 


Now we recall that h\ was the density of the absolutely continuous invariant 
measure of (j)\. But for almost every A € (a, 6), the Bernoulli convolution 
v\ is absolutely continuous and is preserved by 4>\, and so h\ is the density 
of v\. But now any weak* limit point of a family of measures which is ab¬ 
solutely continuous with uniformly bounded density must also be absolutely 
continuous with the same bound on the density. Therefore, since the family 
v\ evolves continuously (in the weak* topology), we see that v\ is absolutely 
continuous for all A G (a, b ) and has bounded density h\. □ 
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4 Computational Techniques 


We have seen in the previous section that showing that 4>\ is piecewise convex 
for all A in an interval would have significant consequences for Bernoulli 
convolutions. Analytically, we have been able to show piecewise convexity 
only for some special values of A for which the distribution F\ is already 
known. We remain optimistic that some further progress could be made here, 
see the comments section. In this section we show how numerical information 
on (J>\ can show convexity up to a certain scale. 

4.1 Showing convexity up to a certain scale for fixed A 

First we choose a natural number M and let X; denote the point jj for 
% G (0,..., M}. We wish to show that 

<P\ {Xi) < ~ (0A (•^'t—l) + (f>\(x i+ 1 )), (9) 

for i < If and A in a certain parameter range. It will then follow that 

<f>\(xj) < ^ Xk ^ + ^ x(yX ^ 

for 0 < i < j < k < This is what we call ‘convexity up to scale jj\ It 
corresponds to the usual definition of convexity restricted to the set of points 
{xo, • • •, xm }, using the fact that 

^ = (Bi) ** + (rii) x ‘- 


Because Lemma 2.2 that tells us that 4>\(x) = 
only need to check ([9]) for i with 1 — A < x* < 


A l x for 0 < x < 1 — A, we 

M 

2 ’ 


For large L we estimate F\(x) by noting that 


Fx(x) < F+ l (x) := 2~ l 


{a\ ■■■ a L G (0, 1} L : (A 1 


L 

1) ajA* < x} 
1=1 
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and 


F x (x) > F x l (x) := 2 L {oi • • • a L G {0, 1} L : (A 1 


= F+(x-\~ l ). 


L 

1) cqA* < x — A l } 
2—1 


Then given the values of F XL (xi ) and F XL (xi ) for z e {1,... ,M}, we can 
bound (p\ from below and above by 

where <p XL and (p XL are defined for 0 < Xi < 1/2 by 

<P\, L ( x i) = V where y is the largest y such that F^ L (y) < 2 F X L (xi), 

< t>x L ( x i ) = V wh ere y is the smallest y such that F XL (y) > 2 F XL {xi). 


Hence, if we have 

< ^((/>x( x i- 1 ) + <t>x (z*+i))> (10) 

for all i with 1 — A<Xj<l/2, then (|9]) holds for all i < Pf. 


The inequalities (10) can be checked with a computer, at least up to errors 
in the floating point arithmetics. We have written a program in CQ that 
calculates the approximations (p XL and (p XL of <j>x,L- Figure [4] shows four 
plots of the approximations, obtained from the mentioned program. 


4.2 Techniques for all A in an interval 


To apply Theorem 3H we would like to show that <px is piecewise convex for 
all A in an interval. We cannot do this computationally, but instead consider 
how to show <px is piecewise convex up to a certain scale for all A in an 
interval. 


We consider a small interval I e = [Ao — e, Ao + e]. The map 

OO 

/ e 9 A4 (A' 1 - 1) ^ 

2=1 

1 The source code for this program is available on the homepage of the second author, 
http: //www. maths, lth.se/matematiklth/personal / tomasp / 
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Figure 4: Plot of 0 AL and 0 AL with L = 24, M = 50, and A = 0.6 (left), 
A = 0.7, A = 0.8 and A = 0.9 (right). 


is differentiable, and if we put D := ( A()+c )(i_ A()+c ) , then 


dA 


(A" 1 -!) J>A* 


i= 1 


- J]a i A i + (A- 1 -l) 

2=1 2=1 

/ _< OO OO \ 

^ i= 1 i \ ' 

/ ^ OO OO \ 

^ i=1 i =1 ' 


< max< — 


< max< — 


= max 


<D, 


. A(1 - A) ’ A(1 — A) I A(l-A) 
holds for all A G I e and all sequences with a % G {0,1}. We conclude that 

L L 

1 - 1) ^ aiA’o < x ==► (A -1 -1)^0^ <x + |A- Ao|E> 


2=1 


2=1 


for any A G I e . Similarly, we have 


(Ao 1 -1) J>a;> 


x 


(A" 1 - 1 )^ 0 ^ >x- |A-A 0 |£ 


2=1 


2=1 


for any A G I e . Using these two estimates we can use F± q l to estimate F± L . 
We get 


f x,l( x ) > f \q,l( x - eD ) and F+ l (x) < F+ Jx + eD). 
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-^0 

e 

convexity to scale 

0.65 

0.000001 

0.02 

0.7 

0 

0.02 

2 -1 / 2 « 0.707106781186548 

0.00001 

0.125 

0.75 

0 

0.02 

0.75 

0.00001 

0.125 

2 -1 / 3 « 0.793700525984100 

0.00001 

0.125 

0.8 

0 

0.02 

0.8 

0.00001 

0.125 

0.85 

0.000001 

0.125 


Table 1: Numerical observations of piecewise convexity to a scale 


Hence, the estimates F t L of F \ 0 gives us estimates on F± L that we can use 
to estimate from below and (p^ L from above. It is then possible to check 
with a computer if the inequalities in ([9]) are satisfied for all A G I e . This has 
been implemented in our program. 

Table [T| shows some result of our program. It displays some values for which 
we have been able to show nummericaly convexity to a certain scale. 

A convolution argument shows that h\ is differentiable for almost all A G 
(2 _ s, 1 ), see [13]. One might suspect that using this information it would be 
possible to show that (p\ is piecewise convex for all A G [2 — 3 , 1 ). However, 
this does not seem to be true, since just as we can sometimes show convexity 
to a scale using numerics, we are sometimes also capable of observing non¬ 
convexity at a certain scale. 

Using our program we have observed that (p\ is not piecewise convex when A 
is the inverse of the root of x 5 + x 4 — x 2 — x — 1 that is larger than 1. We then 
have A ~ 0.8501..., and since 2“3 « 0.7937..., we do not have piecewise 
convexity of (f)\ for all A G [ 2 ~ 3 , 1 ). In this case, 1/A is a Salem number. 

Similarly, the program can be used to show that (j)\ is not convex for A = 
2 1 . Since v\ is known not to be absolutely continuous for this value of A, 
this is not too surprising. We also see a lack of convexity for 1/A equal to 
certain other Pisot numbers. For instance when 1/A is the root of x 4 — x 3 — 1 
or x 3 — x — 1, then (j)\ is not convex to scale 0.005. 
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Let us mention some of the computational difficulties associated with trying 
to prove convexity to a scale for the entire interval I e . Suppose L is even. Our 
program calculates all the 2 L / 2 sums an d stores them in an ordered 

list. This requires quite a lot of memory even for L as small as 60, but the 
time required to preform the calculations is rather short. The sums in the 
list are then combined to get the sums Ylf= 1with double as many terms, 
when needed. This method of storing only the sums of length L /2 instead of 
storing the sums of length L, saves memory but increases computation time. 
However we found that doing so yields a better balance between the use of 
memory and the computation time. 

When A is close to 1, large values of L are needed to get a good accuracy 
in the estimates, requiring an unrealistic amount of memory. This is clearly 
illustrated in Figure |4| where for A = 0.9, the two maps cf>^ L and (p^ L differ 
quite a lot, while for A = 0.6 they are indistinguishable. 

With a computer with 64 GB of memory we are able to run our program for 
L < 60. Table [l] was obtained from running the program with 56 < L < 60. 


5 Further Questions 

There are a number of natural questions that follow on from our work. 

Question 1: Can one show that for all A sufficiently close to 1 we have that 
v\ is absolutely continuous? For almost all A G (2 -1 / n , 1) one has that the 
density h\ is (n — l)-times differentiable, does this extra regularity of the 
density give rise to extra regularity in the functions (j)\! 

Question 2: Can one show that there is an interval J CK containing l/\/2 
such that u\ is absolutely continuous for all A G J. Perhaps this would 
involve showing that the map (j)\ evolves smoothly in a neighbourhood of 
l/\/2. 

Question 3: Are there any other properties of v\ (besides the question of 
absolute continuity) which could be studied using 0 A ? In particular, can one 
forbid the possibility of singular Bernoulli convolutions which have Hausdorff 
dimension 1 by proving a similar result for invariant measures of (p\7 Such 
results do exist in the literature for one-dimensional dynamics, see e.g. [10J 
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and [1], but at present are not in such a form that they would apply to <p\. 
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